A modified Potts model for the interaction of surface-attached polymer complexes 
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We present a simple yet generic model for the behavior of a system of many surface-attached 
flexible polymers with rigid side chains. Beyond its potential application in describing the dynamics 
of the extracellular matrix of mammalian cells, the model itself shows an interesting phase transition 
behavior since the underlying models (a two-dimensional Potts model and a XY-model) undergo 
different phase transitions. 
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I. INTRODUCTION 

Surface-attached polymers are an ubiquitous phe- 
nomenon in living matter. On the cellular level, promi- 
nent examples are the cytoskeletal actin network beneath 
the cell's plasma membrane [l| and the extracellular ma- 
trix 0. While the semiflexible actin polymers interact 
at multiple sites with the intracellular face of the cell's 
membrane, the dominant component of the extracellu- 
lar matrix, the hyaluronic acid (HA), rather is a flexi- 
ble polymer, end-grafted to the extracellular face of the 
membrane and extending to the cell's environment. To 
provide an efficient barrier that is capable of protecting 
the cell, extracellular HA is modified by rigid aggrecan 
combs that are linked to the flexible HA backbone. As 
a consequence, the cell is protected by a 'jungle' of HA- 
aggrecan complexes. 

We have previously shown that attaching a single rigid 
side chain to a flexible, end-grafted polymer (resembling 
the HA-aggrecan complex) can lead to a considerable 
stiffening of the backbone [3j , hinting at a mechanism by 
which cells can tune the rigidity of their protective ma- 
trix. In this study, we had used a dynamic Monte Carlo 
algorithm on a regular cubic lattice to study the steady- 
state properties of a flexible, end-grafted backbone poly- 
mer with a rigid side chain. Extending this approach 
to a system of many interacting backbones, however, is 
computationally challenging due to the massive increase 
in the system's autocorrelation time. 

An alternative approach to a brute-force Monte Carlo 
sampling is a further coarse graining of the system that 
reduces the problem to its basic observables, namely the 
height of the backbones and the orientation of the rigid 
side chains. Here, a combination of the Potts model (de- 
scribing the height levels) and the two-dimensional XY- 
model (describing the orientation of the side chains) ap- 
pears as a promising candidate. Using this approach, 
however, one has to deal with an intricate problem of 
statistical physics - the prediction of the system's phase 
behavior and phase transition. The Potts model (with 
q < 4 height levels) undergoes a second order phase tran- 
sition [4J while the two-dimensional XY-model shows a 
Kosterlitz-Thouless transition 0, @. It is far from ev- 



ident which phase transition will be encountered when 
coupling these two systems. Aiming at modeling the be- 
havior of the extracellular matrix in terms of such a com- 
bined model thus requires to elucidate the model's phase 
behavior in the first place. 

Inspired by the potential value of the combined Potts- 
XY-model in describing features of the extracellular ma- 
trix, we have investigated the phase transition of the 
model. In particular, we have studied the phase transi- 
tion in a model that couples the q = 3 Potts model with 
the XY-model by means of extensive Monte Carlo simu- 
lations. We find that the phase transition is dominated 
by the Potts model, i.e. a second order phase transition is 
observed. The ordered state at low temperatures shows 
domains of parallel aligned XY-spins with the Potts lev- 
els collectively exploring all states in a stochastic fashion. 
In contrast, for high temperatures the XY-spins are dis- 
ordered with an average Potts level (q) = 1 that only 
varies to a minor extent. Relating these findings to the 
original biological problem, we speculate that changing 
the interaction strength between individual HA-aggrecan 
complexes, i.e. driving the system through the phase 
transition, may be one way for the cell to build up a 
stochastically fluctuating protective barrier. 



II. MODEL 

As described previously Q , the basic unit of the extra- 
cellular matrix may be reduced to a single flexible, end- 
grafted backbone with a rigid side chain (Fig.[T]). While 
a simulation of this unit yields valuable insights into the 
stiffening of the individual backbone due to the attached 
rigid side chain, many of the details may not be necessary 
to understand the generic behavior of the multi-complex 
system. The basic quantities describing the interaction of 
many individual units are the height p of each side chain 
above the substrate and its orientation angle 9 with re- 
spect to the y-axis. Thus, the individual backbone side- 
chain complexes are interpreted as mutually interacting 
rotors aligned on a regular two-dimensional lattice. 

We assume here that the sites at which the backbones 
are attached to the substrate are sufficiently far aparat 
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FIG. 1: For describing as system of end-grafted, self-avoiding 
polymer backbones with rigid side chains (left), only the es- 
sential degrees of freedom, i.e. the height p of the side chain 
above the substrate and its orientation angle 9, are retained. 
A system of mutually interacting complexes is obtained by 
using a two-dimensional lattice of such rotors (right) where 
the dynamics is determined by Eq. {3). 



from each other so that any interaction between the com- 
plexes is mediated by the side chains only. This also 
means that side chains only interact when they have the 
same height p. This aspect can be modeled by a Potts- 
type Hamiltonian 0, 

H P = -K P J2 S PiPi ! fte{0,l,...,g-l}. (1) 

(id) 

Here, pi denotes the height of side chain i above the 
substrate, q is the total number of allowed heights, and 
Kp > is the coupling constant (in units of (3 — fc^T). 

The interaction between side chains with the same p 
will be considered as a nearest-neighbor ferromagnetic 
interaction leading to a prefered parallel alignment of the 
side chains. This kind of interaction seems plausible as a 
first approach since it mimics the repulsive forces between 
side chains that may be due to steric or short-ranged 
electrostatic potentials. Consequently, we use for this 
part of the interaction the two-dimensional XY-model 
Hamiltonian: 

Uxy = -K X y S * S i = - RxY s < s i cos % > ( 2 ) 

{■id) (id) 

with the coupling constant Kxy > and denoting 
a summation over nearest neighbors. For simplicity, the 
length of the rotors (=side chains) Si — \Si\ is set to 
unity while the angle Oij denotes the relative orientation 
of the side chains i and j. 

Combining the above models Eq. (fl]) and Eq. ^ yields 
the hybrid model that is used to describe the surface- 
attached polymer complexes of the extracellular matrix: 

H = -K ^2 Spipj ^iSj cos 9ij ; pi e {0, . . . , q - 1} . 

(id) 

The two-dimensional Potts model is known to undergo 
a second order phase transition for q < 4 |4J while the 
XY-model shows a Kosterlitz-Thouless transition H Q . 
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FIG. 2: The inner energy U per rotor shows a sigmoidal in- 
crease with temperature T for the pure XY-model [Eq. J2J], 
the pure 3-state Potts model [Eq. ([TJ] , and the hybrid model 
[Eq. Q] . The transition temperature T c for the latter appears 
to be shifted towards lower values while the slope in the tran- 
sition region is strongly increased. Full lines are guides to the 
eye. 



For q = 3 our model Eq. ^ can be interpreted as a 
3-state Potts model with a randomly varying coupling 
strength Kij = cos Oij between neighbouring Potts spins 
i and j. As this local coupling may become very small 
and even zero, the nature of the phase transition of the 
hybrid model may be influenced by vacancies that are 
known to drive the transition of the (q < 4)-state Potts 
model to a first order behavior Q. Given this variety of 
phase behaviors, it appears difficult to predict the sys- 
tem's behavior a priori. 



III. SIMULATION METHOD 

For simulations of the hybrid model Eq. we relied 
on the Monte Carlo (MC) method (see, e.g., [lfj for an 
introduction). Far away from the critical temperature, 
we employed the Metropolis MC algorithm In each 
step, a rotation of each rotor (side chain) by some angle 
A9i and a change of its Potts level (height) were proposed 
and this combined move was accepted or rejected for each 
rotor according to the Metropolis criterion. 

Close to the critical temperature, the standard 
Metropolis MC suffers from critical slowing down, i.e. 
it becomes difficult to generate enough statistically inde- 
pendent configurations for a reliable statistics. To over- 
come this problem, we used a cluster MC approach [l2j], 
i.e. a modification of Wolff's cluster algorithm [l3|. Here, 
one deals with two degrees of freedom, one discrete (Potts 
level p) and one continuous (orientation angle 6). Thus, 
two different kinds of clusters are grown and flipped in 
each step - one with respect to p (embedded in an en- 
semble of fixed Oi's) and another one with respect to 8 
(embedded in an ensemble of fixed pi's). 

The treatment of the Potts degree of freedom has 
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to be carried out carefully as ferromagnetic and anti- 
ferromagnetic bonds can occur between neighboring ro- 
tors. In the beginning, a pair of values for p is chosen 
randomly, i.e. (0,1), (0,2), or (1,2), respectively. Ro- 
tors in the remaining Potts state are kept unchanged, 
i.e., the identity operator is applied, and thus do not con- 
tribute to the ongoing cluster growing step. The effec- 
tive (local) coupling between rotors i and j is K e s = 
K ■ SiSj cos 9i j. For K e g > a ferromagnetic bond is es- 
tablished with probability P — 1 — exp (— K e g5 Pi ^ Pj ) be- 
tween rotors on the same Potts level. If K e g < 0, an anti- 
ferromagnetic bond is formed with probability-? = 1 — 
exp (K e g (1 — 8 PuPj )) between rotors on different Potts 
levels. When the cluster is grown, all Potts levels of 
the involved rotors are switched according to the initially 
chosen pair of p- values (e.g. «-» 1). 

Following Wolff's embedding trick for growing and flip- 
ping clusters in the XY-model [l3| the rotors are pro- 
jected onto the :r-axis. This results in an Ising model 
of the re-components Sf with random couplings. These 
'effective' spins are connected by a bond with probability 
P = 1 - exp {5 PuP] ■ min[0, 2KSfSj}) . Please note that 
the Potts states of the rotors are accounted for by S PiiP . 
meaning that bonds are only formed between rotors on 
the same Potts level. Finally, when the cluster is grown, 
all ^-components of the constituing rotors are inverted. 



IV. RESULTS 

To study the thermodynamic properties of a system of 
interacting rotors modelled by the Hamiltonian Eq. ([3|) 
with ? = 3we conducted extensive MC simulations using 
the method described in the preceeding section. The lin- 
ear system size L (i.e. using N = L 2 rotors) was varied 
in the range 32 < L < 1024 in order to apply a finite-size 
scaling analysis. We have concentrated here in particu- 
lar on the inner (potential) energy U per rotor and the 
transition temperature T c . 

Metropolis MC simulations of small systems revealed 
significant differences in the behavior between the hy- 
brid model and the two constituting models (Potts and 
XY). Fig. [2] shows the variation of U with temperature 
for all three models. Obviously, the transition point T c 
of the hybrid model is much lower than those of the con- 
stituting models. Furthermore, the transition region is 
much steeper, suggesting a first-order transition. At this 
point, however, it cannot be determined if the slope really 
diverges at T c which would be a criterion for a discon- 
tinuous phase transition. This behavior would be qual- 
itatively different from that shown by the constituting 
models, yet would agree with the random Potts model 
Q . Please note that for T < T c the energy U of the hy- 
brid model resembles essentially that of the correspond- 
ing pure XY-model. Thus, the Potts degree of freedom 
seems to be 'frozen' in this regime while all excitations 
occur with respect to the XY-degree of freedom. 

For a more thorough investigation, we used the cluster 




FIG. 3: (a) The inner energy U below and above the critical 
temperature T c shows only a single maximum for a system 
with N = 256, while for T ~ T c a slight double-peak structure 
is observed, (b) Indeed, this double-peak structure at T ~ T c 
is even better visible for small systems and seems to subside 
for larger system sizes. 



algorithm described in the preceeding section. In con- 
trast to the standard Metropolis MC method it allowed 
for extensive simulations of large systems close to T c due 
to a reduced dynamic exponent. We focussed on the inner 
energy U as we found no proper order parameter to follow 
the phase behavior of the hybrid model. As can be seen 
from Fig. [3^, the distribution P(U) shows pronounced 
single peaks for T«r c and T » T c , while for T « T c a 
double-peak seems to emerge. In Fig. Eb, the distribution 
P{U) is shown for three different system sizes at T s=s T c . 
The distributions were obtained by reweighting [3] of 
distributions close to T c . Our numerical data indeed 
highlight a double-peak structure in P(U), indicating the 
coexistence of two phases and thus a first-order phase 
transition. To extrapolate if this finding persists in the 
thermodynamic limit, we first considered the logarithm 
of the relative gap depth A = ln[P(U) max /P(U) m i n ] at 
T w T c for growing system sizes L. For a first-order 
transition, this quantity is supposed to grow with L. 
In the hybrid model studied here, however, A — > for 
L — > oo (Fig. Hh.). Thus, the energy distribution assumes 
a Gauss-like shape in the thermodynamic limit, consis- 
tent with a second-order phase transition with only a sin- 
gle phase existing in the system at T c . Considering the 
finite-size scaling of the transition temperature T c with 
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FIG. 4: (a) The depth A = ln[P(U) m ax/P(U) min } of the 
minimum between the double peak in P(U) appears to van- 
ish with system size L, thus indicating a second-order phase 
transition, (b) In agreement with this, the phase transition 
temperature shows a scaling T c ~ 1/L. (c) The tendency of 
the 4* -order cumulant towards a value V = 2/3 is also in 
support of a second-order phase transition. 



the inverse system size 1/L supports this notion since we 
find T c oc 1/L as anticipated for a second-order transition 

(Fig. Hp). 

Further evidence for a second order phase transition 
comes from the 4 th order cumulant of the energy distri- 
bution [H 
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FIG. 5: Snapshot of the system at (a) T < T c and (b) T > 
T c . Rotor orientations are indicated by arrows, Potts levels 
p — 1,2,3 are labeled in black, dark grey and light grey, 
respectively, (c) Probability distribution of the mean Potts 
level in the system P((p)). In the ordered state at T < T c 
(dashed), all Potts states p — 1, 2, 3 are populated with similar 
probabilities while for T at T c (dash-dotted) and T > T c (full) 
the state p = 1 is clearly favored. 



structure are characterized by V ^ §. While a small 
dip in V(T) is observed for small system sizes (Fig. 0b), 
it subsides for large L and only very small deviations 
from V = 2/3 are seen. In the thermodynamic limit, one 
thus expects V = I, i.e. a Gauss(-like) distribution for 
all temperatures. 

The phase behavior of the system above and below the 
critical temperature is shown in Fig. [5] While for T < T c 
the rotors (=rigid side chains) are ordered (Fig. they 
assume random orientations for T > T c . The mean Potts 
level (=thc average height of the rigid side chains) in 
both cases is (q) = 1 (Fig. [5b), yet for T < T c collective 
fluctuations of the rotors to the three degenerate Potts 
levels are observed. Given the small transition probabil- 
ities for these height fluctuations for fixed temperatures 
(which we circumvented via the MC cluster algorithm) 
the dwell time at an individual Potts levels may be very 
large. 



V. CONCLUSIONS 



V = 1 



3- {u 2 y 



(4) 



For a Gaussian distribution, this quantity takes on the 
value V — | while distributions with a double-peak 



In this work, we have proposed a model for a system 
of surface-attached polymers with rigid side chains that 
is found, e.g., in the extracellular matrix of mammalian 
cells. To reduce the computational effort when simulat- 
ing large systems of end-grafted flexible polymers with 
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rigid side chains, a coarse grained approach has been 
chosen combining two well-known models of statistical 
physics, namely the Potts and the XY-model. The ther- 
modynamics of the 'hybrid' model, given by the Hamilto- 
nian Eq. ([3]) , has been studied in Monte Carlo simulations 
and a second-order phase transition was found. 

Given our results, it seems reasonable to call the hybrid 
model rather a 'modified Potts model' than a 'modified 
XY-model' as the phase behavior of the system is domi- 
nated by the 3-state Potts model: At low temperatures, 
the rotors essentially condensate into a single Potts state. 
This corresponds to the (three-fold degenerate) ground 
state of the Potts model showing long-range order with 
respect to p. Excitations in this regime arise due to per- 
turbations in the ordered state of the XY-model. At high 
temperatures, the influence of the XY degree of freedom 
becomes manifest as a randomized coupling of the Potts 
spins which includes the existence of vacancies that may 
drive the transition of the 3-state Potts model from sec- 
ond to first order 0. Indeed, at first glance, the tran- 
sition of the hybrid model appears to change character 
(cf. Fig. [5] and Fig. [3]). Only more thorough investiga- 
tions revealed the transition to be still continuous, i.e. to 
remain second order. Thus, the vacancies evoked by the 
random coupling drive the transition towards first order, 
but do not suffice to fully change its character. 



Coming back to the original problem that inspired the 
model, i.e. the dynamics of the extracellular matrix, the 
observed ordering of the system for T <T C may actually 
be used by a cell to change the thickness of its protective 
layer. A cell may use the transition to the regime T <T C 
to freeze the extracellular matrix in a desired homoge- 
nous ordered state/height that is not necessarily 'thicker' 
than in the disordered regime. Thus, speaking in terms 
of the original problem, the cell may tune the thickness 
of its extracellular matrix by tuning the interactions of 
the (polyelectrolytic) aggrecan side chains, thereby driv- 
ing the system through the phase transition. It will be 
interesting to experimentally manipulate the aggrecan in- 
teraction by changing ion conditions for the extracellular 
matrix and thereby to induce the above described phase 
transition. 
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